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ABSTRACT 

We present the discovery of substellar-mass companions to three stars by the 
ongoing Penn State - Toruh Planet Search (PTPS) conducted with the 9.2-m 
Hobby-Eberly Telescope. The K2-dwarf, BD +14 4559, has a 1.5 Mj companion 
with the orbital period of 269 days and shows a non-linear, long-term radial 
velocity trend, which indicates a possible presence of another planet-mass body 
in the system. The K3-giant, HD 240210, exhibits radial velocity variations that 
require modeling with multiple orbits, but the available data are not yet sufficient 
to do it unambiguously. A tentative, one-planet model calls for a 6.9 Mj planet 
in a 502-day orbit around the star. The most massive of the three stars, the K2- 
giant, BD +20 2457, whose estimated mass is 2.8+1.5 M , has two companions 
with the respective minimum masses of 21.4 Mj and 12.5 Mj and orbital periods 
of 380 and 622 days. Depending on the unknown inclinations of the orbits, the 
currently very uncertain mass of the star, and the dynamical properties of the 
system, it may represent the first detection of two brown dwarf-mass companions 
orbiting a giant. The existence of such objects will have consequences for the 
interpretation of the so-called brown dwarf desert known to exist in the case of 
solar-mass stars. 

Subject headings: planetary systems-stars: individual (HD 24210, BD +14 4559); 
brown dwarfs: individual (BD +20 2457) 
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Introduction 



Searches for planets around giant stars offer an efficient way to extend studies of plan- 
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the most massive sub stellar companions to early-type stars are possible (IGalland et al.ll2005 



Lagrange et al.ll2009l ). it is much more efficient to utilize the power of the radial velocity 
(RV) method by exploiting the fact the GK-giants, the descendants of the main sequence 
A-F type stars, have cool atmospheres and sufficiently many narrow spectral lines to make 
achieving a <10 m s _1 RV measurement precision possible. 

The GK-giant surveys are beginning to provide the statistics, which are needed to 
constrain the efficiency of planet formatio n as a function of ste llar mass and chemical com- 
position. In fact, the initial analyses by iJohnson et al.l (120071 ) and lLovis &: Mayor! (120071 ) 
extend to giants the correlation between planetary masses and the masses of their primaries 
observed for the lower-mass stars. Most likely, this is simply because massive stars tend 
to have more massive disks. These results are in accord with the core accretion scenario 



of planet formation ([Kennedy fc Kenyonl 120081 ). Furthermore, iPasquini et al.l (120081 ) have 
used the apparent lack of correlation between the frequency of planets around giants and 
stellar metallicity to argue that this effect may imply a pollutio n origin of the observe d 
planet frequency - metallicity correlation for main sequence stars (IFischer &: Valentil 120051 ) . 
Similar conclusions h ave been reached in a preliminary study of our giant star sample by 
dzielihski et~al1l2009h . 



The steadily extending baselines of the ongoing surveys of post-MS giants begin fur- 
nishing multiplanet system detections that are needed to study the dynamical evolution of 
planets around the off main-sequence stars. The multiplic i ty of planetary sys tems around 
MS stars has been firmly established (jUdry fc Santos! 120071 ; IWright et al.ll2009l ) and there is 
no reason why this tendency should not be present among giants. In fact, the ongoing search 
for planets around GK-giants with the Hobby-Eberly Telescope (HET) by our group (Penn 
State- Toruh Planet Search, hereafter PTP S) has already tentativ ely identified a two-planet 
system around the KO-giant HD 102272 (jNiedzielski et al.ll2009l). an d long-term trends in 



the RV data have been seen in other surveys (e.g. ISato et al.l (j2008al )) 



In this paper, we describe the detection of a planet around the K3-giant, HD 240210, 
and show evidence that the star has more low-mass companions, and the intriguing discovery 
of two brown dwarf-mass bodies orbiting the K2-giant, BD +20 2457. We also describe the 
discovery of a Jupiter- mass planet and a non-linear RV trend in the K2-dwarf, BD +14 4559. 
This detection is the first result of the extension of the PTPS to evolved dwarfs in the upper 
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envelope of the MS. This part of the project will be described in detail in a future paper. 

The plan of this paper is as follows. An outline of the observing procedure and a 
description of the basic properties of the three stars are given in Section 2, followed by the 
analysis of radial velocity measurements in Section 3. The accompanying analysis of rotation 
and stellar activity indicators is given in Section 4. Finally, our results are summarized and 
further discussed in Section 5. 



2. Observations and properties of the stars 



Observations were made with the Hobby-Eberly Tele scope (HE T) (IRamsey et al.lll998l ) 
equipp ed with the High Res olution Spectrograph (HRS) (lTulllll998l ) in the queue scheduled 
mode (IShetrone et al.l 120071 ). The spectrograph was used in the R=60,000 resolution mode 
with a gas cell (I2) inserted into the optical path, and it was fed with a 2 arcsec fiber. 
Details of our s urvey, the observing procedure, and data analysis ha ve been described in 
detail elsewhere (INiedzielski et al.l 120071 ; iNiedzielski fc Wolszczanl 120081 ) 



Radial velocities were measured using the standard I2 cell calibration technique (IButler et al. 



19961 ) . A template spectrum was constructed from a high- resolution Fourier Transform Spec- 
trometer (FTS) I2 spectrum and a high signal-to-noise stellar spectrum measured without 
the I2 cell. Doppler shifts were derived from least-squares fits of template spectra to stellar 
spectra with the imprinted I2 absorption lines. The radial velocity for each epoch was de- 
rived as a mean value of the independent measurements from the 17 usable echelle orders 
with a typical, intrinsic un certainty of 6-8 m s -1 at lcr-level. This RV precision levels made 
it quite sufficient to use the IStumpfil (119801 ) algorithm to refer the measured RVs to the Solar 
System barycenter. 

The long-term precision of our RV measurements has been verified by the analysis of 
data derived from the monitoring of stars that do not exhibit detectable RV variations. The 
results for the K0 giant BD+70 1068 ( V=9.29, B-V=1.15 ± 0.04, T e// = 4580 ± 17 K) are 
shown in Figure 1 as an example. The relative RVs for this star vary with a=12 m s -1 , 
whereas our estimated intrinsic RV uncertainty for this series of observations was 7 m s _1 . 
The excess RV scatter is consistent with the ap proximate 10 m s" 1 amplitu de of solar-type 
oscillations derived from the scaling relations of iKjeldsen fc Bedding! (119951 ). 



The atmospheric parameters of the stars under consideration were taken from Zielinski 
et al. (2009, in p r eparat ion, hereafter Z09), who have estimated their values using the method 
of Takeda et al. (|2005aB l. With these values and the luminosities estimated from available 
data, stellar masses were derived by fitting the ensemble of parameters characterizing the 
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star_(log(L/L ), log T e //, log (g), and metallicity) to the evolutionary tracks of lGirardi et al. 
( 120001 ) . The parameters of the three stars are summarized in Table 1. 



2.1. BD +14 4559 

Radial velocities of BD +14 4559 were measured at 43 epochs over the period of 1265 
days between MJD 53546 and 54811. Typically, the signal-to-noise ratio per resolution 
element in the spectra was 150-260 at 594 nm in 10-25 minutes of integration, depending on 
the atmospheric conditions. The estimated mean RV uncertainty for this star was 8 m s _1 

BD +14 4559 (AG +14 2370, HIP 1 04780 LTT 16221) is a high p roper motion K5 
rtTuron et all (|l993h ) sta r with V=9 m 63 JWeisI ll986h . B-V=0 m .98±0.08 JrlQg et all lioooh 
and 7r= 19.99+1.61 mas (jPerryman fc ESAl 119971 ). T he trigonom etric parallax is consistent 
with the photometric value of tt= 24 mas derived by (jWeislll986l ). The atmospheric parame- 
ters of BD +14 4559 as determined in Z09 indicate that it is a dwarf with log(g)=4.60±0.05, 
T e// =5008±20 K, and [Fe/H] =0.10+0.07. 



Using the existing Hipparcos (jPerryman fc ESA 1997) and 2MASS ( Skrutskie et al. 



20061 ) photometry and the empirical calibrations of (jRarmrez fc Melendezl 120051 ). we have 
estimated the effective temperature of the star to be T e y/=4814±26 K. As this value of T e ff 
is too high for a K5 dwarf, we conclude that the spectral type of BD +14 4559 is K2V. 

The derived absolute bolometric magnitude of the star, My, and log(L/L ), are 5.56 
and -0.32, re s pectiv ely. By comparing the stellar parameters to evolutionary tracks of 
Girardi et all J2000h . we have computed the mass of BD +14 4559 to be M/M =O.86 ± 
0.15 and R/R©= 0.95 ± 0.2. 

The projected rotational velocit y of BD +14 4559 , vs ini =2.5 ± 1 km s _1 , was estimated 
using the cross-correlation method (IBenz fc Mayor! Il984l ). From this value and the adopted 
stellar radius we have obtained an estimate of the rotation period of P rot =19 days, which is 
much shorter than the observed 267-day period of the RV variations. From the uncertainty 
of the vsin(i) determination alone, the stellar rotation period may range from 13 to 32 days. 



2.2. HD 240210 



Radial velocities of HD 240210 were measured at 38 epochs over the period of 1655 days 
between MJD 53187 and 54842. Typically, the signal-to-noise ratio per resolution element 
in the spectra was 150-300 at 594 nm in 4-10 minutes of integration, depending on the 
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atmospheric conditions. The estimated mean RV uncertainty for this star was 8ms 



V= 



HD 240210 (BD +56 2959) i s a K7 (ICaimon fc Pickerinelll924f) star with V=8 m 33, B- 
l m .The parallax of 6.5+0.04 (fade et alJ |200oh and tt= 7.6+8.7 mas JPerrvman fc ESA 



19971 ) are consistent with the 7r= 7.0±2.6 mas of (IGatewood et al.lll992l ). The atmospheric 
parameters of HD 240210 (Z09) indicate that the star is a giant with log(g)=2.31±0.11, 
T e// =4297±25 K, and [Fe/H]=-0.18±0.12. 

1993) and 2MASS fekrutskie et aL 



Using the existing Hipparcos (jPerryman fc ESA 



20061 ) photometry, and the empirical calibrations of (jRamfrez fc Melended l2005h . we have 
estimated the effective temperature of the star to be Tg//=4290±13 K. These values of 
effective temperature and log(g) are consistent with a K3-giant star. 

The derived absolute bolometric magnitude of the star and logfL/L^), are 0.38 and 



1.75, respectively. By comparing stellar parameters to evolutionary tracks of iGirardi et al. 
f)2000h we have computed the mass of HD 240210 to be M/M =1.25 ± 0.25 and R/R©= 13 
± 3. 

The projected rotational velocity of HD 240210 , vsini < 1.0 ± 1 km s _1 , was estimated 
using the cross-correlation method. From this value and the adopted stellar radius we have 
obtained an estimate of the rotation period of P rot > 654 days. 



2.3. BD +20 2457 



Radial velocities of BD +20 2457 were measured at 37 epochs over the period of 1833 
days between MJD 53033 and 54866. Typically, the signal-to-noise ratio per resolution 
element in the spectra was 140-260 at 594 nm in 12-30 minutes of integration, depending on 
the atmospheric conditions. The estimated mean RV uncertainty for this star was 7 m s _1 



BD +20 245 7 (AG +20 1166) is a K2 flHeckmannlll975f ) star with V=9 m 75, B-V=l m .25±0.09 
(Ho g et al.l 120001) w hich is similar to V=9.57 and B-V=1.35 as measured photometrically by 



( jBusko et al.lll980l ). who observed it as a comp arison star in their stell ar variability study. 
The parallax of the star, 7r= 5.0+26.50 mas (jPerryman fc ESAl Il997t ). is obviously very 



poorly determined. The atmospheric parameters of BD +14 4559 (Z09) indicate that it is a 
dwarf with log(g) = 1.51±0.05, T e// =4137±10 K, and [Fe/H]=-1. 00+0.07. 



Using the existing Hipparcos (jPerryman fc ESA 1997) and 2MASS ( Skrutskie et al. 



20061 ) photometry and the empirical calibrations of (jRamfrez &: Melendezl 120051 ). we have 
estimated the effective temperature of the star to be T e jj=4127±17 K. We conclude that 
the spectral type of BD +20 2457 is K2II. 
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The derived absolute bolometric magnitude of the star, My, and log(L/L ), are - 
3.0 and 3.17, respec tively. By comparing the stellar parameters to evolutionary tracks of 



Girardi et al.l (120001 ). we have estimated the mass of BD +20 2457 to fall in the 1.3 - 4.3 
M/M range (2.8±1.5 M/M ), and R/R©= 49 ± 26. 

As for the other two stars, the projected rotational velocity of BD +20 2457 , vsini <1.0 
km s _1 , was estimated using the cross-correlation method. From this value and the adopted 
stellar radius we have obtained an estimate of the rotation period of the star, P rot > 2460 
days. 



3. Analysis of the radial velocity data 



The process of modeling the RV measurements of the three stars discussed in this 
paper is documented in Figures 2-6. In the first step of the analysis, single, 6-parameter, 
Keplerian orbits w ere least-squares fitted to data using the Levenberg-Marquardt algorithm 
JPress et al.lll992h . For the K-giants, HD 240210 and BD +20 2457, the estimated, 7-8 m s" 1 
errors in RV measurements were clearly insufficient to account for the observed RV variations 
in the post-fit residuals. This excess RV variability is believed to originate i n fluctuations of 



the s tellar surface and radial and non-radial oscillations of the giant stars (jde Ridder et al. 



20091 ) with the typical RV amplitudes around 20 m s 1 (IHekker et al.ll2006l ). To account for 
this effect, we have quadratically added 10 m s _1 and 30 m s -1 to the calculated RV errors 
for HD 240210 and BD +20 2457, respectively. 

Evidently, in all the three cases, the best-fit residuals from single orbit models leave non- 
random trends that need to be modeled further. A statistical significance of these trends 
was assessed by calculating false alarm probabilities (FAP) for a null hypothesis that the 
trends can be adequately accounted for by a single Keplerian orbit and noise. Th e FAPs were 
calculated with the aid of the radial velocity scrambling method (jWright et al.l (120071 ). and 
references therein). As illustrated in Figure 6, FAP < 0.1% for the three data sets, which 
indicates that the additional, systematic RV variations seen in the data are statistically 
significant. 



3.1. BD +14 4559 

As shown in Figure 3, the best-fit, single Keplerian orbit model for the RV variations of 
this star leaves behind a significant, non-linear trend, which reaches the amplitude of ~50 m 
s _1 over the 1300-day span of observations. At this point, rather than attempting to search 
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the x 2 ' space for a range of acceptable two-planet solutions, we have modeled the observed 
RV variations with a Keplerian orbit and a parabolic trend added to it. The best-fit x 2 value 
for this model drops from 4.4 for the single orbit case to 2.1 indicating, together with the 
FAP <0.1%, that the trend is significant. 

The final solution (Table 1) includes a planet with the minimum mass of 1.5 Mj up , in 
a 269-day, 0.78 AU, e =0.29 orbit around the star. Given the stellar mass of 0.86 M Q , the 
mass function of 2.5 x 10~ 8 M , and the above RV amplitude and time span of the observed 
trend, one can broadly constrain a putative second companion to be another planet with the 
minimum mass >2.4 Mj up and the orbital radius >2.3 AU, assuming e =0. 

The remaining post-fit rms residual, ctrv= 11 m s _1 , slightly exceeds the estimated 8 
m s _1 precision of our RV measurements for this star. This excess "jitter" is within the 
estimat ed, intrinsical ly generated RV noise of 5 m s _1 due to solar-type oscillations in stable 



dwarfs (IWrighti 120051 ) . 



3.2. HD 240210 

A single orbit fit to the RV data for this star (Figure 4, Table 2) gives x 2= 9.8, the post- 
fit rms residual, cr^y=39m s _1 , and the very clear unmodeled RV variations. Our search 
for a multiple orbit solution has resulted in several two- and three-planet models, which 
gave nearly identical improvements of the fit. This is illustrated in Figure 4 by a two-planet 
model, which produces x 2 =5.0, (Jrv=25 m s -1 , FAP < 0.1%, and the persistent lower-level 
systematics in the RV residuals. Because the existing measurements are obviously insufficient 
to obtain satisfactory constraints on the Keplerian model of the RV variations of this star, 
we list in Table 2 the initial one-planet solution and postpone further analysis until the time, 
when enough data become avaiable. The provisional, best-fit orbit involves a 6.9 Mj up planet 
in a 502-day orbit with the semi-major axis of 1.3 AU and e =0.15. 

We also note that the 25 m s _1 post-fit residual from the multiplanet fits is about 3 
times larger than the formal precision of the RV measurements for HD 240210. As mentioned 
above, this is in agreement with the estimated RV scatter in K-giants due to their internal 
activity. 



3.3. BD +20 2457 

A single, Keplerian orbit fit to the RV data for this star is shown in Figure 4. The resid- 
uals from this fit are characterized by % 2 =13.7, £Trv=105 hi s _1 , and a ~600-day periodicity 
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with the amplitude in excess of 100 m s 1 . Once again, as shown in Figure 6, the FAP for 
the existence of the second periodicity is less than 0.1%. 

As shown in Figure 5, the best-fit, two-orbit solution offers a dramatic improvement 
to the Keplerian model of the RV variations observed in BD +20 2457, with the values of 
X 2 and (Try reduced to 5.2 and 60 m s -1 , respectively. The model parameters calculated 
for the stellar mass of 2.8 M and listed in Table 3 call for a system of two substellar-mass 
companions with the respective minimum masses of 21.4 Mj up and 12.5 Mj up , orbital periods 
of 380 and 622 days, semi-major axes of 1.45 and 2.01 AU, and mild eccentricities of 0.15 
and 0.18. 

As in the case of HD 240210, the estimated, 7 m s _1 uncertainty of the RV measurements 
of BD +20 2457 is much smaller than the actual, 60 m s" 1 rms of the post-fit residuals. 
This additional RV "jitter" broadly agree s with the 138 m s~ l ampli tude of the solar-type 



oscillations extrapolated for this star from lKjeldsen fc Bedding! (119951 ). 



4. Stellar photometry, rotation, and line bisector analysis 



In order to verify that the observed RV periodicities are indeed caused by the Keplerian 
motion, we have thoroughly examined the existing photometry data in search for any periodic 
light variations, and performed a complete analysis of line bisectors and curvatures for each 
of the three s tars discussed in this paper. This analysis used the cross-correlation method 
proposed by ([Martinez Fiorenzano 1 120051 1 . For each star, the cross-correlation functions 



were computed from ~1000 line profiles with the I? lines rem oved from the spectra. This 
procedure has been described in detail in iNowak et al.l (120091 ). 



Because both the Ca II K emission line and the infrared Ca II triplet lines at 849.8-854.2 
nm are outside the range of our spectra, we have used the Ha line as a chromospheric activity 
indicator. As the shape of the Ha line in the spectra that showed significant telluric line 
contamination could be affected by it, care was taken to omit such spectra from the analysis. 
Equivalent widths (EW) of the line profile, defined as I/I c <0.9, were measured within the 
central 1 A of the absorption profile, where the signature of stellar activity should be most 
pronounced (the uncertainty of the EW defined this way depends on the actual depth of the 
Ha line and the signal-to-noise ratio). 
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4.1. BD +14 4559 



The existing photometric databases for this star cover a wide range of epochs from, 
MJD 47891 to 54792, and include ex tensive time series of measureme nts from the Hipparcos 
and T vcho JPerrvman fc ESAl ll997h. the NSVS JWozniak et alil2004h . and the ASAS exper- 



iment ( Pojmanski 19971 ). None of the four data sets reveal any periodic light variations. In 



particular, the ASAS data, which are contemporaneous with our RV measurements of BD 
+ 14 4559, include observations made at 280 epochs between MJD 52754 and MJD 54792. In 
this case, the light variations are characterized by the mean value of V=9.650±0.021 mag, 
which is consistent with the earlier measurements. As an example, the Lomb-Scargle (LS) 
periodogram of the ASAS photometry data in shown in Figure 7. 

We have obtained 43 bisector and line curvature measurements, which resulted in the 
mean bisector velocity span, BVS=4.6 ± 32 m s -1 , and the mean bisector curvature, BC=-8 
± 27 m s _1 with no trace of any periodic variability. No correlation between the BVS, BC, 
and RV was found (r=0.09 and r=0.08, respectively). Similarly, our analysis of 23 spectra of 
the star has shown that the EW of the H Q line exhibited a 3% rms scatter around its mean 
value of 396 ± 12mA, which was not correlated with the observed RV variability (r = -0.15, 
while the critical value of the Pearson correlation coefficient at the confidence level of 0.01 
is r 21i0 .oi=0.53). The LS periodograms of the H a EW , the BVS, and BC for BD +14 4559 
are shown in Figure 7. Clearly, the only observable that exhibits periodic time variations in 
time is the radial velocity. 

Using the scatter seen in the ASAS photometry of the star, and its rotational velocity 
determined above, we can estimate the amplitude of radial v elocity variat ions and of the 
BVS due to a possible presence of a spot on the stellar surface (IHatzes M2002I ) . The observed 
radial velocity amplitud e of BP +14 4 559 is three times larger than the 35 m s _1 total RV 
amplitude predicted by IHatzes I (120021 ) for a spot with a filling factor of f=0.02, on a star 
rotating at vsini =2.5 km s _1 . Within the estimated uncertainty of the rotation velocity, 
the expected RV amplitude induced by a hypothetical spot is at least two times less than the 
observed one. The expected bisector variations of 5 m s" 1 are comparable to the precision of 
our RV measurements and cannot have a detectable effect on our results. To summarize, the 
facts that there is no variability in the available photometric measurements over the period 
of ~19 years, the observed period of the RV v ariations is much longer than the estimated 
rotation period of the star, a significant eccentricity must be included in the best-fit model 
of the RV curve, and that the observed RV variations have been consistent over the five 
consecutive cycles of the period, make their interpretation in terms of a rotating spot very 
unlikely. 
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4.2. HD 240210 

In the case of this star, the o nly long-term photometr y available comes from 129 Tycho 



measurements with the Vr filter (jPerryman &: ESAj 119971 ) . between MJD 47946 and 49016 



These data give the mean Vr=8.499 ± 0.012 and the scatter in Vj- of 0.136 mag. No peri- 
odic variability has been detected in these measurements. In ad dition, there are 19 ep ochs 



of photometric observations of this star available from the NSVS (jWozniak et al.1 120041 ). col- 
lected between MJD 51362 and 51493. These data are characterized by the median V=7.965 
± 0.01 and the scatter of 0.038. Since the NSVS photometry is of a much better quality 
than the Tycho data, we consider this scatter more believable. As shown in Figure 8, the 
LS periodogram of the Tycho photometry of HD 240210 reveals no significant periodicities. 

The mean values of the BVS and the BC for this star are 62 ± 23 m s _1 , and 23 ± 
27 m s, respectively, with no statistically significant variability. No correlations between the 
RV residuals and the BVS and BC were found. The respective correlation coefficients are 
r=0.08 for the BVS and r=-0.18 for the BC (r 36) o.oi=0.40). Furthermore, the average EW of 
the H a line measured over 28 epochs amounts to 602 ± 13, and the EW itself does not show 
any periodicities within 2%. We have considered a possible rotational modulation in H a 
by calculating the Pearson correlation coefficient with the RV data residuals after removing 
the putative orbit. We found no correlation between the H a EW and the RV residuals (r= 
-0.26, whereas r 2 6,o.oi = 0-48). The LS periodograms of the four observables for HD 240210 
discussed above are shown in Figure 8. As in the case of BD +14 4559, radial velocity is the 
only one that varies periodical ly. 

The observed photometric scatter in HD 24210, if interpreted in terms of a spot rotating 
with the star would also be detectable in both the radial velocity and the line bisector data. 
The estimated total RV variations should be of the order of 23 m s _1 and the bisector 



variations should amount to 31 m s 1 (IHatzes II2002I ). The semi-amplitude of the observed 



periodic signal in radial velocity is much larger than the predicted variations induced by 
a hypothetical spot of size consistent with the existing photometeric variations. Periodic 
variations in line bisectors produced by a hypothetical spot should manifest themselves as 
the uncertainty in our BVS measurements, which is comparable the the predicted value. 
In any case, these variations should show a periodicity equal to the stellar rotation period, 
P > 654 d. Because no such variations are present in our bisector measurements for this 
star, we conclude that the observed scatter is a combination of low quality photome try and 
possible solar type oscillations of the star. Such oscillations might account for 3.3 mmag in 
the photometric variability and V osc = 13 m s _1 in the BVS using the scaling relations of 



Kieldsen & Beddind (jl995f l. 
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4.3. BD +20 2457 

As in the case of BD +14 4559, we have identified four photometric data sets that we 
could utilize in ou r analysis, but we have concentrated on the data from the ASAS project 



(jPojmanskil 119971 ) . which contains extensive measurements that are contemporaneous with 
our RV observations. These data include 120 measurements made between MJD 52622 and 
MJD 53400 and the 191 additional ones made between MJD 52623 and MJD 54911. We have 
used 319 grade A measurements from the ASAS to derive the mean V=9.693 ± 0.051 mag 
and to compute their LS periodogram in search for any periodicities. As shown in Figure 
9, no significant periodicities are present in the spectrum. Similarly, no periodicities were 
found in the other three, lower quality time series. 

We have obtained 37 bisector and curvature measurements for this star, with the mean 
BVS=105 ± 39 m s -1 , and the mean BC=47 ± 30 m s -1 , both revealing no significant 
periodicities in the corresponding LS periodograms (Figure 9). We have also searched for 
correlations between the RV residuals and the BVS and BC after having removed the putative 
companions b and c. The respective correlation coefficients are r=-0.33 and r=-0.19 for 
BVS and r=-0.12 and r=-0.09 for BC, respectively (the critical value of Pearson correlation 
coefficient at the confidence level of 0.01 is r35 i o.oi=0.42). The EQ of the H Q line for BD 
+20 2457 is 653+ 23, and this parameter shows no significant periodicities within 3.5% for 
23 observations (Figure 9). As for the other two stars, we have searched for a possible 
rotational modulation in H a by calculating the Pearson correlations coefficient with the RV 
data residuals after removing the individual orbits. We found that there is no significant 
correlation between the H Q EW and the RV residuals (r=-0.21 and -0.51, respectively, and 
r 2i,o.oi=0.53). 

A photometric scatter resulting from a spot rotating with BD +20 2457 would induce 
the estimated total RV variations on the order of 30 m s" 1 and the bisector variations of 40 



m s _1 (IHatzes Il2002l ). The semi-amplitude of the observed periodic signal in radial velocities 
is over one order of magnitude larger. As in HD 240210, the anticipated variability in 
line bisectors produced by a periodic spot rotation is similar to that observed in our BVS 
measurements, these variations should show a periodicity equal to the star's rotation period, 
P > 2640 d, which appears not to be the case (our observations cover only 1833 days of the P 
> 2640 day rotation period). Since no measurable variations are present in the bisectors and 
line curvatures, we conclude that the observed photometric scatter is again a combination 
of low quality photometry and possible solar type oscillations of the star. Such oscillations 
might account for 5 . 9 mm ag and up to V osc = 138 m s _1 using the scaling relations of 



Kieldsen & Beddind (J1995J). 
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Discussion 



The most recent results from the Penn State- Toruh search for extrasolar planets de- 
scribed in this paper add a Jupiter-mass planet to the large body of planets around main- 
sequence stars (lSantosll2008l) and increas e the s ti ll modest number of substellar-mass object s 



around GK-giants ( e.g. iDollinger et al.l (120091) : iLiu et all (120091 ): iNiedzielski et all (J2009J); 



Sato et al.l (j2008bl ); INiedzielski et al.l (120071 ) and references therein) by at least three new 
detections. 

In principle, the long period RV variations in red giants may be related to non-Keplerian 
effects, including the stellar r otation modulation of surface inhomogeneities, intrinsic activity, 
or non-radial pulsations (e.g. lHatzes et al.ll2006l ). Our analysis of the available photometric 



data and the behavior o f line bisectors in the three stars, following the established prac- 
tices (jOueloz et al.ll200ll ). has shown no significant correlations between the time variability 
of these stellar activity indicators and the RV variations. Consequently, our RV measure- 
ments find the most plausible explanation in terms of the Keplerian motion of substellar- 
mass companions arou nd the observed stars. The details of this analysis can be found in 
Niedzielski et all toO^ ). 



In the case of the K2-dwarf, BD +14 4559, the RV data reveal a planet with the minimum 
mass of 1.5 Mj in a 269-day, 0.78 AU, e=0.29 orbit around the star. An additional, long-term 
trend seen in the data suggests the presence of another, long-period companion, which is very 
likely to be another planet. This adds to the growing number o f confirmed and anticipated 
multiplanet systems around solar- type stars (jWright et al.ll2009l ). 



A single planet model of the RV variations in the K3-giant, HD 24210, leaves highly 
correlated post-fit residuals signaling a possible presence of additional substellar-mass bodies 
orbiting the star. The presently available data are not sufficient to obtain an unambiguous 
multibody solution for this system. The provisional parameters for one planet that can be 
fitted for give a 6.9 Mj body in a 502-day, 1.33 AU, e=0.15 orbit that will have to be revised, 
when another planet (or planets) are added to the current model. Nevertheless, the solution 
is robust enough to conclude that HD 24210 has at least one planet with large minimum 
mass. This result further emph asizes the observed c orrelation between th e stellar mass and 
the masses of orbiting planets (I Johnson et al.l 120071 : lLovis &: Mayorl 120071 ). 



The RV variations observed in the K2-giant, BD +20 2457, can be unambiguously 
modeled by two Keplerian orbits with the periods of 380 and 622 days, 1.4 and 2 AU semi- 
major axes, and the respective eccentricities of 0.15 and 0.18. For the estimated stellar mass 
of 2.8 M , these parameters yield minimum masses of the orbiting bodies of 12.5 Mj and 
21.4 Mj, respectively. The ~ 1.6 ratio of the orbital periods places the orbits relatively close 
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to the 3:2 mean motion resonance (MMR). These characteristics suggest that the two bodies 
share a common origin from a disk that once surrounded the primary star. At the same time, 
at least one of its companions, and very possibly both of them, have masses that exceed the 
13 Mj deuterium burning limit that conventionally separates planets from brown dwarfs. 

Further analysis of this intriguing system is complicated by the fact that, due to the 
highly uncertain parallax of the primary, its mass is poorly known and spans a range that 
is at least as wide as 1.3-4.3 M Q (Section 2). This is illustrated in Figure 10, which shows 
the constraints on the masses of the two substellar companions under the usual assumption 
of a random distribution of orbital inclinations. There is a 95% probability that the masses 
of the two bodies fall in the respective ranges delimited by the inclinations of 18° and 90°. 
Clearly, high masses, approaching the hydrogen burning limit are very unlikely. Similarly, 
a possibility that both masses stay below the deuterium burning limit is marginal at best, 
unless one accepts an unrealistically low primary mass, approaching 1 M or even less. 
Consequently, within the currently available constraints, it is reasonable to assume that the 
true masses of the two companions to BD +20 24 57 have masses in the brown dwarf range. 
An extensive study of the dynamics of this system will be described in a subsequent paper. 

The two brown dwarf - mass bodies discovered by our survey around BD +20 2457 add 
to the previous three detections of single objects of this kin d around intermediate-mass stars 
(IHatzes et al.l 120051 ; lLovis fc Mayorl 120071 ; iLiu et al.l 120081 ). Only one such tw o-companion 



system has been previously identified around a solar-mass star, HD 168443, by iMarcy et al. 
(120011 ). On the basis of the facts presented above, it appears entirely reasonable to assume 



that the two companions to this star have originated from a massive disk and acquired enough 
mass to place them above the formal 13 Mj brown dwarf limit. This raises an intriguing 
possibility that, in t he case of substellar-m ass companions to giants, both the concept of the 
brown dwarf desert (IMarcy fc Butlerll2000l ) and the meaning of the deuteriu m burning limit 
will h ave to be revisited. In fact, this problem has been remarked on by lLovis fc Mayor 
(120071 ) in the context of the observed correlation between the stellar and the planetary mass. 
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Table I. Stellar parameters of BD + 14 4559, HD 24210 and BD +20 2457. 



Parameter 


BD +14 4559 


HD 24210 


BD +20 2457 


V 


9.63 


8.33 


9.75 


B-V 


0.98+0.08 


1.65+0.04 


1.25+0.09 


Spectral type 


K2V 


K3III 


K2II 


7r [mas] 


19.99+1.61 


7.0+2.6 


5.0+26 


Te// [K] 


5008+20 


4290+13 


4127+17 


log g 


4.60+0.05 


2.31+0.11 


1.51+0.05 


[Fe/H] 


0.10+0.07 


-0.18+0.12 


-1.00+0.07 


log L/L 


-0.32 


1.75 


3.17 


M*/M 


0.86+0.15 


1.25+0.25 


1.3-4.3 



Table 2: Orbital parameters of the BD +14 4559 planet. 



Parameter 


BD +14 4559 b 


P [days] 


268.94 ± 0.99 


T [MJD] 


53293.71 ± 6.35 


K [ms- 1 ] 


55.21 ± 2.29 


e 


0.29 ± 0.03 


lu [deg] 


87.64 ± 7.87 


rri2 sin % [Mj] 


1.47 


a [AU] 


0.777 


(xl) 1/2 


2.10 


a RV [ms x ] 


11.43 
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Table 3: Tentative orbital parameters of the new planet around HD 240210. 



Parameter 


HD 240210 b 


P [days] 


501.75 ± 2.33 


T [MJD] 


54486.76 ± 10.93 


K [ins" 1 ] 


161.89 ± 3.49 


e 


0.15± 0.02 


uj [deg] 


277.49 ± 7.77 


m2 sin % [Mj] 


6.90 


a [AU] 


1.33 


(xl) 1/2 


9.79 


a RV [ms" 1 ] 


38.9 



Table 4: Orbital parameters of the sub-stellar-mass companions to BD +20 2457 derived 
from the best fit of a two-body Keplerian model to the RV data. 



Parameter BD +20 2457 b BD +20 2457 c 



P [days] 


379.63 ± 2.01 


621.99 ± 10.20 


T [MJD] 


54677.03 ± 28. 


19 53866.95 ± 27.99 


K [ms" 1 ] 


322.35 ± 9.57 


160.03 ± 7.21 


e 


0.15 ± 0.03 


0.18 ± 0.06 


uj [deg] 


207.64 ± 21.99 


126.02 ± 16.54 


m 2 sin i [Mj] 


21.42 


12.47 


a [AU] 


1.45 


2.01 


(xl) 1/2 




5.20 


a RV [ms -1 ] 




60.02 
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54000 



Fig. 1. — Relative radial velocities for KO-giant BD+70 1068. The star shows RV variability 
at a— 12 m s" 1 whereas the intrinsic RV measurement uncertainty for these observations was 
7 ms -1 . The remaining RV scatter can be attributed to solar - type oscillations at the level 
of 10 m s _1 
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Fig. 2. — Top: Radial velocity measurements of BD +14 4559 (filled circles), the best-fit of a 
single planet Keplerian model to data (solid line) and the fit of that model with a parabolic 
trend (dashed line). Center: The post-fit residuals for a single planet model. Bottom: As 
above for the planet with the trend fitted out. 
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Fig. 3. — Top: Radial velocity measurements of HD 240210 (filled circles), the best-fit single 
planet Keplerian model (solid line), and a tentative, best-fit model of two planets to data 
(dashed line). Center and bottom: See Figure 2. 
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Fig. 4. — Top: Radial velocity measurements of BD +20 2457 (filled circles) and the best-fit 
of a single planet Keplerian model to data (solid line). Bottom: The post-fit residuals for 
the above model. 
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Fig. 5. — A decomposition of the observed RV variations in BD +20 2457 (filled circles) 
into contributions from the orbital motions of two brown dwarf-mass companions. From top 
to bottom: The observed radial velocities (filled circles) and the best-fit Keplerian orbit for 
companion b (solid line) with the contribution from companion c fitted out. The same for 
companion c with companion b fitted out. The best fit of both b and c to the RV data. The 
post-fit residuals from the two-companion fit. 
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Fig. 6. — Histograms of the values of a/x^ obtained from the fits of the Keplerian models 
to scrambled RVs used to estimate the FAPs for the three stars discussed in the text. In 
each case, 1000 sets of scrambled RVs have been generated. Vertical arrows point to the 
respective \J~xl values derived from fits to unscrambled velocities. In all three cases the FAP 
values are less than 0.1%. 
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Fig. 7. — The Lomb-Scargle periodogram s of (a) Radial ve locities (b) Bisector Velocity Span 
(c) H Q equivalent width and (d) AS AS (jPojmanskil 1 199 71 ) V photometry of BD +14 4559. 
The levels of FAP=1.0 % and 0.1% are shown. 
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Fig. 8. — The Lomb-Scargle periodograms of (a) Radial velocities (b ) Bisector Velocity 
Span (c) H a equivalent width and (d) Tycho (jPerryman &: ESAlll997l ) Vt photometry of 
HD 240210. The levels of FAP=1.0 % and 0.1% are shown. 
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Fig. 9. — The Lomb-Scargle periodogram s of (a) Radial ve locities (b) Bisector Velocity Span 
(c) Hq, equivalent width and (d) ASAS JPoimanskilll997h V photometry of BD +20 2457. 
The levels of FAP=1.0 % and 0.1% are shown. 
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Fig. 10. — Constraints on the masses of substellar companions to BD +20 2457. The solid 
lines delimit the masses of the inner companion computed from the mass function over the 
estimated 1.3-4.3 M range for the stellar mass and orbital inclinations from 18° to 90°. The 
dashed lines have the same meaning for the outer companion. The dotted lines constrain 
the theoretical mass range for brown dwarfs. 



